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Abstract 



We present a new class of hermitian one-matrix models originated in the 
Wi+00 algebra: more precisely, the polynomials defining the W"i+oo gen- 
erators in their fermionic bilinear form are shown to expand the orthogonal 
basis of a class of random hermitian matrix models. The corresponding po- 
tentials are given, and the thermodynamic limit interpreted in terms of a 
simple plasma picture. The new matrix models can be successfully applied 
to the full bosonization of interesting one-dimensional systems, including all 
the perturbative orders in the inverse size of the system. As a simple appli- 
cation, we present the all-order bosonization of the free fermionic field on the 
one-dimensional lattice. 
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The study of one-dimensional fermionic systems has recently attracted much at- 
tention due to progress in manufacturing of mesoscopic devices. Although much 
progress has been made in the theoretical understanding of the dynamics of these 
strongly correlated systems, many interesting questions remain still unanswered. 
One theoretical tool that is successfully used is the bosonization of the fermionic 
degrees of freedom [|I|: for a given fermionic system in the thermodynamic limit, 
one considers the low-lying particle-hole excitations around the Fermi points. Lin- 
earization of the dispersion relation around these leads to a relativistic conformal 
field theory of these gapless effective degrees of freedom. More recently, an ex- 
tended bosonization procedure incorporating the subsequent non-linear corrections 
of the actual dispersion relation has been proposed 0, ||. The algebraic bosoniza- 
tion technique is an all-order bosonization approach based on the Wi+oo algebra, a 
characteristic fermionic symmetry extending the conformal symmetry of the stan- 
dard bosonization. The main advantage of this method is that it yields a systematic 
expansion of the full hamiltonian (including interactions) and all physical observ- 
ables in powers of the momentum fluctuations around the Fermi points, involving 
algebraic procedures only. 

In this letter, we present a one-parameter family of hermitian one-matrix models 
which originate in the Wi+oo algebra and are intimately connected to the algebraic 
bosonization of one- dimensional fermionic systems. Consider the Fock space of a 
(1 + l)-dimensional relativistic Weyl fermionic field with creation and annihilation 
operators al^ and a„, respectively (n is an integer that labels the momentum), sat- 
isfying {al, am} = ^n,m- Wc define the VTi+oo operators (A; = 0, 1, 2, . . .) in this 
specific basis 0: 



[ Vi ] = {jn - zm)K+^-i + q{t, j, n, m)V^^' + ■■■ + 6^^n+^mfl c rf(z, n) , (2) 



where the structure constants q{i, j,n,m) and d{i,n) are polynomial in their ar- 



k = 0, 1, . . . , [(« + j — l)/2], where the brackets denote the integer part function 
(the complete expression of (H) is given in 0). The meaning of the operators Vj^ 
is as follows: they parametrize particle-hole excitations near the Fermi points, such 
that the index n accounts for the fixed momentum transfer, whereas the index k 
denotes the geometric moments (i.e.,"multipole moments") of the deformation of 
the Fermi-Dirac sea 0, |^, ^. The explicit expression of the polynomials PJ^ {x) can 
be inferred directly from the Wi+oo algebra 




(1) 



which satisfy the Wi^, 



■oo 



algebra 0, 
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. ) -° ^ ^ 

where (a)„ = n"=o^('^ + ^) denotes the Pochhammer symbol. Some properties can 
be deduced from (^): i) the polynomials have a definite parity with respect to the 
reflection about x = 1/2: even n (resp. odd) correspond to even (odd) polynomials. 
a) for n large enough, the zeroes of -P^(x) are located on the line x = 1/2 + with 
C real, and for pure imaginary values of n, all zeroes fall along this line. For k even, 
all zeroes come in complex conjugate pairs, and for k odd the same is true except 
for a single zero at x = 1/2. However, the most important consequence of is 
that the Pn{x) satisfy a three-term recurrence relation: 

^P^ix) = Pr^(x) + 1 P^{x) + Pt\x) . (4) 

This property has the important consequence of leading to a new class of random 
matrix models associated to the Wi+oo algebra, as we shall show in the following. 

Matrix models constitute another important tool when studying systems dis- 
playing universal behavior in the thermodynamic limit [^. There are many know 
examples in which the symmetries that characterize the relevant effective degrees of 
freedom of a given system are simply encoded in a specific random matrix model, 
which posses the advantage of being a simpler (0 + l)-dimensional theory. Consider 
a generic random one-matrix model, defined by its partition function: 

Z{N) = J d^' ^ exp [-TrV(<l>)] , (5) 

where $ is an x Hermitian matrix, Tr is the trace, ane the real function V{^) 
is the potential. Eventually, we shall be interested in the limit N oo, which is the 
relevant thermodynamic limit. The matrix $ has real eigenvalues Xi, i = 1, . . . , N, 
in terms of which the partition function can be written as 0: 

Z{N) = n{N) j . . . j l[dx, A^\x) exp(^-J2y{x,)j , (6) 

where An{x) = Yli<ji^j - ^i) and n{N) = (27r)^(^-i)/V nf=i ^' is the volume of 
the "angular" (non-diagonal) entries of $ P]. A general strategy to compute Z{N) 
is to factorize the A^ integrals on the r.h.s. of (^), which can be done provided there 
exists a set of orthogonal polynomials Vk{x) satisfying: 

/oo 
dx e-^(") Pk{x) Pi{x) = hk 5k,i . (7) 
-oo 
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In turn, this condition can be met provided the polynomials Vk{x) satisfy a three- 
term recurrence relation given by: 

xVk{x) = Vk+i{x) + Sk Vk{x) + Rk Vk-i{x) , (8) 

where Sk and Rk are specific to the matrix model A crucial equality that relates 
the potential with the recursion relation is Rk = hk/hk-i, and the partition function 
is given by: 

N-l N~l 

Z{N) = Q{N) N\ l[hk = Q{N) N\ J] . (9) 

k=0 k=0 



In the following, we shall focus in the set of polynomials (^, and show that they 
are orthogonal with respect to an specific measure which is then used to define the 
associated matrix model. One anticipates, therefore, an infinite family of models 
parametrized by the momentum n, and defined by the family of potentials Vn{x). 
Given the location of the zeroes of the polynomials (^), it is natural to perform a 
"Wick rotation" in x and define V!^{x) = i^P^ {—ix + (n + l)/2) as the set of VTi+oo 
orthogonal polynomials. They posses a definite parity under refiection in x; the first 
few are: Vl{x) = 1, Vl{x) = x, Vl{x) = x^ + {n^ - 1)/ 12 ,V^{x) = x^ + x{3n'^ -7)/20. 
The recursion relation (^) now takes the form: 

xV^ix) = V^^\x) + ^l^f^-Pt'ix) , (10) 

which implies the existence of the set of potentials Vn{x) making them orthogonal. 
These functions are readily found once one recognizes that the polynomials V^ix) 
are Hahn polynomials, for which the orthogonality relations are known [p!0|. There 
are, in fact, two classes of Hahn polynomials: discrete and continuous. The first 
case arises for n integer, such that there is a finite set of polynomials, since the 
recursion relation ( p!OD stops ai k = n. However, given (P, we are interested in the 
case for which the family is infinite: this is the continuous case, that arises for n 
not an integer. For our purposes, it is convenient to consider the values n = ip, 
p = 0,1,2,... (i.e., imaginary momentum) since this yields a positive and non- 
degenerate inner product. This condition can be viewed as the Wick rotation in 
momentum space that corresponds to the one previously performed in x. Under 
these conditions, the orthogonality relations are given by [|T^: 

r dx g,{x) Vl{x) VUx) = \nk + l + ^p)\^^ 



k 

TT 

9p{x) = TTT^ TT T • (11) 

cosh(/7ra;j -|- cosh(p7rj 
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Figure 1: The potentials Vp{ function of x for p = 0, 1, 2. Note the asymptotic 

hnear behavior discussed in the text. 

From this expression, one extracts the family of potentials Vp{x) defining the matrix 
models as in fBI): 



- In qJx) 



In 



cosh(27ra;) + cosh(p7r) 



TT 



p = 0,1,2,... , (12) 



which asymptotically behave as Vp{x) ^ 27r|2;| — 21n(27r) for \x\ > 7p/7r, jp = 
ln[27rcosh(p7r)] (see figure 1). From (pT]), one also knows that the norm of the k — th 
polynomial in the p ~ th matrix model is: 



ht 



\T{k + l + tp)\' 



i2k + l] 



2k 



2 ' 



(13) 



which is in agreement with the recursion relation (p!0| ) through the equality Rp = 
hp/hp~^ = k'^{k'^ + p^)/4(4fc^ — 1). From ([13D , one can therefore define the partition 
function Zp{N) and free energy density fp{N) for the p — th model as: 



N-l 



Zp{N) = n{N) iV! n ^5 = exp{N'fp{N)) . 

fc=0 

The asymptotic form of fp{N) in the thermodynamic limit ^ oo is: 

fp{N) 

independent of p to the given order. 



(14) 



^ In AT - I'^ + iln^) + O (In AT/AT) 



(15) 
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The thermodynamic hmit leads also to a simple physical picture in terms of a 
fictitious one-dimensional plasma: consider the partition function for the poten- 
tial ([I^ ) (we consider for simplicity the case p = 0). The asymptotic form of Vq{x) 
is linear, and suggests the definition of a new scaling variable y = x/N, appropriate 
in the limit N oo. This implies that Vo{x) ^ 2TTN\y\ — 21n(27r) for \y\ > 1/A, 
with A = ttN/'-jq. In the scaling limit, y is defined in the entire real axis, except 
for a small region of the size of the cutoff 1/A around the origin. We can write the 
partition function approximately as: 



N 



I n dy, exp -nUttY^M-I^Y.^ Hv^ " Vjf + f{N) 

(16) 

where f{N) = (A^ - 1) In + 2 ln(27r). In the large A^ limit, Zq{N) is dominated by 
the saddle-point configurations {yi, . . . ,yN} satisfying the equations: 

2vr sg(y,) = I- 5^ , t = l,...,N, (17) 

where sg{y) is the sign function. The scaled variable y is appropriate for taking 
the continuum limit: yi — >• y, ^.(. . .)/N — > J dyp{y){. . .), where p{y) is the density 
function. The saddle-point equations are: 

vr sg(2/) = V ^ dy' , (18) 

J-oo y-y 

p{y')dy' = 1, (19) 

where V denotes the Cauchy principal part of the integral, and the second equality is 
the normalization condition. These equations give a semiclassical physical interpre- 
tation to our system: they describe a system of fictitious equally charged "particles" , 
located on a one-dimensional rod, and subject to a neutralizing background. The 
particles repel each other with a two-dimensional Coulomb force, while they are at- 
tracted with constant force to an infinite two-dimensional uniformly charged plane, 
that is perpendicular to the rod. One can find a solution to equations ([T9| ) by as- 
suming that the system is symmetric under y — ^ —y, and therefore the total number 
of "particles" is even, with A^/2 of them on each side of the infinite charged plane. 
Therefore, one seeks for a solution of the form p{y) = p+{y)E{y) + p_{y)E{—y), 
with E{y) the Heaviside function. In the following we shall focus on one side only, 
e.g. (+). As it is customary, one introduces the function 

F4y) = r dy' , (20) 

J-oo y-y 

defined on the complex half-plane Re {y) > 0. Assuming that this function has 
one cut running along the real axis from to a, one can solve (|T9|) by a standard 
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procedure yielding: 



= vr -n J^-^ , (21) 



From (pID one determines the density: 



pM = (22) 

where a = l/vr is determined from the normahzation condition (|1^). One can verify 
that along the cut and close to its extreme ?/ ~ a, the behavior of p+{y) is in 
agreement with the expected one for a linear background potential |p. 

Let us now focus on the application of the previous results to one-dimensional 
condensed matter models. Consider a free non-relativistic fermionic field on a one- 
dimensional lattice of unit spacing. We shall outline the effective field theory for- 
mulation of this problem in terms of the VTi+oo symmetry generators While 
this description does not appear as necessary in the free theory, it becomes rather 
advantageous when including the effects of the interactions, as, for example, in the 
XX Z Heisenberg model. We shall postpone the full discussion of the effects of the 
interactions for a forthcoming publication, but remark here that the scope of our 
approach is to develop a bosonization scheme valid to all orders in the expansion pa- 
rameter, which is the inverse size of the system, given by N. In the thermodynamic 
limit, the non-relativistic field can be expanded in the basis of the relativistic Weyl 
field, and the Fermi sea interpreted as a Dirac sea [Q. Moreover, the free hamilto- 
nian H can be written as the sum of two independent terms, H+ and describing 
the fluctuations around the right (-I-) and left (— ) Fermi points at momenta in^?, 
respectively, with np = (iV — 2)/4; i.e., we consider the ground state of the system to 
be at half-filling, for simphcity. Both terms are decoupled, up to global (zero-mode) 
constraints. Therefore, we shall consider only in the following: 

n+ = - ^ f{nF + r) -.alar. , f{n) = cosi—n] , (23) 

r=— 00 ^ ^ 

where are the fermionic Fock operators of the Weyl field, obtained from the 
original ones by shifting the momentum index, and r is integer. The dispersion 
relation around the Fermi point given by f{np+r) is non-linear, and can be expanded 
in terms of the zero- momentum W^i+oo polynomials (^ as [Q: 

00 

7f+ = 5^ Cfc V,' , (24) 
where the coefficients are such that: 



sm 



2n / 1 

— r 

V 2 



E ^'^ ^oir) . (25) 



k=0 
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Making the change of variables r —>■ {—ir + 1/2) and passing to the basis VqIv) 
allows us to determine using the orthogonality relation (|^): 

Cfc = y dx go{x) sin \^—rj V^{x) A; = 0, 1, 2, . . . , (26) 

where for convenience we have made the redefinition 27[/N ^ i2TT/N, so as to deal 
with regular expressions; this analytic continuation should be undoed at the end of 
the computations. From ([26| ) one immediately concludes that C2n = for n integer, 
and defines Ck = C2k+i {k = 0,1,2,.. .) for the other cases, which can be worked out 
by considering the set of integrals: 

/ + 00 
dxg,{x) e'^yV^ix) , n = 0,1,2, ... . (27) 
-oo 

As a consequence of (p!0D, they satisfy the recursion relations: 

i,Uv) - -^'-^ - K i.-dv) . 'M-^-i^y (28) 

which are formally solved by /„(?/) = Vq (^—i-^^ hiv)- The coefficients (^6]) are 
then readily determined yielding: 

(-1)'' /27t\ 
Ck = 7^ W(^^J , A: = 0,l,2,.... (29) 

As and example, after analytically continuing back rr/N — > —m/N as explained 
above, 



sin(7r/iV) 



1 - 



A^tan(7r/iV) 




which reproduces the results of 0. We have also verified the leading term in C2 = 
— (27r/A^)^/6 + . . . yields the correct result. 

Finally, we remark that the matrix models considered in this letter belong to the 



class of models with non-polynomial potentials, which include also the Penner ||12 
and Frustrated Spherical |T^ models. According to the 't Hooft correspondence, one 
can associate to each matrix model a tiling of a two-dimensional random surface, 
such that the the number of sides and the ways of gluing the elementary poligons are 
determined by the algebraic powers of the monomials that define the potential [|ll[] 0. 



The class of matrix models with non-polynomial potentials possess the characteristic 
property of having an unbounded set of polygons and gluing vertices among them. 
We conclude by remarking that it would be interesting to explore whether other 
infinite-dimensional Lie algebras ecould give rise to further new classes of random 
matrix models. 

GZ acknowledges the hospitality of the Abdus Salam Center for Theoretical 
Physics (Italy) and the Physics Department and I. N.F.N, of Florence (Italy). 
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